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Wc investigate convergence of the density matrix renormalization group (DMRG) in the thermo- 
dynamic limit for gapless systems. Although the DMRG correlations always decay exponentially in 
the thermodynamic limit, the correlation length at the DMRG fixed-point scales as £ ~ m 1 ' 3 , where 
m is the number of kept states, indicating the existence of algebraic order for the exact system. The 
single-particle excitation spectrum is calculated, using a Bloch-wave ansatz, and we prove that the 
Bloch-wave ansatz leads to the symmetry E(k) — E(tt — k) for translationally invariant half-integer 
spin- systems with local interactions. Finally, we provide a method to compute overlaps between 
ground states obtained from different DMRG calculations. 



I. INTRODUCTION 

Since White constructed the density matrix renor- 
malization group (DMRG) technique about five years 
agca'B, numerical renormalization techniques have be- 
come very useful. The DMRG has by now been applied 
to a wide range of different problems beyond quantum 
spin systems, where it was originally used. Today, peo- 
ple use it to compute, for instance, -properties of two- 
dimensional classical lattice systemstl thermodynamics 
of one- dimensional quantum ^stemso etc. The DMRG 
and matrix product stateaj E9 have proven to be com- 
putationally very efficient and to determine properties of 
many systems with unusually high accuracy. 

This paper aims at a better understanding of the un- 
derlying structure and the fundamental limitations of the 
DMRG. It has been reported that the DMRG is les 
curate for gapless systems than for gapped systemstiJ'l 
This has motivated us to analyse the DMRG of gapless 
systems in more detail, which we do in two steps: first, 
we investigate the correlation functions (which pertain to 
the wavefunction) , and second, we study the excitation 
spectrum. We have choosen to study free fermions on a 
one-dimensional lattice as a paradigm of gapless systems. 

In this work, we show that the DMRG of a gap- 
less system converges, for each choice of the number of 
kept states, to a fixed-point. The corresponding corre- 
lation functions are calculated by using an eigenvalue 
techmqueEj'Lj, and the results show that this fixed- 
point describes with increasing accuracy the exact sys- 
tem. In particular, we address the question of how the 
DMRG handles algebraic correlations (infinite correla- 
tion lengths), which is characteristic of gapless systems, 
in contrast to gapped systems having finite corrrelation 
lengths and determine a scaling formula for how the 
DMRG correlation length depends on the number of kept 
states. In DMRG applications, this scaling formula may 
serve as a guide for how many states that is neccessary 



to keep in order to accuratly calculate correlations. We 
perform calculations for the particle-hole and the density- 
density correlation functions. In addition, we derive con- 
ditions for which types of operators that can give truly 
long ranged DMRG correlations. By introducing a gap, 
we also investigate how the DMRG correlation lengths 
change as the critical point (free fermions) is approached. 

A matrix product Bloch-wave ansatz .-has been pro- 
posed for describing the excited statesEJ'EJ of a system. 
In this paper, we have used the Bloch-wave ansatz to 
calculate the excitation spectrum. In particular, we 
look at the spectrum close to the Fermi points, where 
the gap closes. Furthermore, it was recently shownlij 
that many half-integer spin systems have the symmetry 
E{k) = E(ir — k). We prove that this symmetry is inher- 
ent also in the Bloch-wave ansatz. 

In DMRG calculations, it is important to check con- 
vergence with respect to the only source of error (ex- 
cept from round off errors), namely the truncation of the 
Hilbert space. A commonly used measure of the trti 
cation error is the truncation of the density matrixo 
However, that measure is algorithm dependent and there- 
fore not universal. It would be desirable to instead cal- 
culate and use the overlap between states from differ- 
ent DMRG calculations as a measure. The problem is 
then that the states will refer to differently renormal- 
ized Hilbert spaces. In this work, we demonstrate how 
the matrix product formalism can be used to handle this 
problem. 

The organization of the paper is as fol lows : the system 
we have studied is defined in Section II A ; a brief in- 



troduction to matrix prod uct s tates and the Bloch-wave 
ansatz is given in Section II B; we decribe how we cal- 
cula te co rrelation functions and overlaps in Section II C 
and II D; computational methods are outlined in Section 

f. The results are presented and discussed in Section 
convergence of the projection operator (to a fixed- 
point) and of the ground state is demonstrated in Section 
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Ill A ; correlation functi ons and the scaling formula are 
treated in Section [II B; Appendix |A| contains an impor- 
tant result for the determination of correlation lengths in 
a fermionic system; conditions for true long range corre- 
lations are derived in Appendix |bJ the excitation spectra 
is presented in Section III C and a proof of the symmetry 
of the spectra is given in Appendix Final ly, t he results 
and conclusions are summerized in Section IV 



II. METHODOLOGY 



A. Hamiltonian 



We have studied a system of non-interacting spinless 
fcrmions on a one-dimensional lattice. The Hamiltonian 



is: 



N 



h.c.] 



j=i 



N 

t 



(i) 



where N is the size of the lattice, cj creates a fermion on 
site j, and t is the hopping amplitudcES (t = 2 throughout 
this work). We have added a staggered on-site potential 
e to the lattice since we want to compare DMRG of the 
gapless system (e = 0) to that of gapped systems (e ^ 0). 
The simple Hamiltonian gives us the advantage of having 
access to exact solutions when evaluating the DMRG. In 
the remaining part of this section, we will state exact 
results used in our analysis. 

The Hamiltonian is particle-number conserving and is 
invariant under the transformation Cj — > (— l) J c| +1 for all 
e. In addition, when e = 0, the Hamiltonian has particle- 
hole symmetry, i.e. it is invariant under the transforma- 
tion 

We will only consider chains of length TV = An + 2 in 
order to have a unique-ground state, which corresponds 
to a half filled systemO. The gap between the valence 
band and the conduction band at the Fermi points is 
2\e\. The correlation functions C(l) decay algebraically 
for the gapless system and exponentially for a gapped 
system. When e = 0we have 



C ph {l) = (c]c i+i ) = -^sin7r//2, 



and 
C M (l) 



ynj+i) - (nj}(n j+ i) 



(2) 



sin 2 ttI/2 (3) 



for the particle-hole and density-density correlation func- 
tions respectively, where rij = CjCj. 

When e ^ 0, the correlation length for the ground 
state particLe.-hole correlation function can be calculated 
analyticallyEJ to be 



In e/t + y/1 + {e/tf 



(4) 



Similarly, the exact density-density correlation length is 
given by £ dd = £ ph /2. 

Finally, we note that there is a well known connection 
to spin-systems. We associate a spin-1/2 with each site 
in the lattice and consider an occupied site as spin up 
and an empty site as spin down. Using the spin raising 
and lowering operators we may rewrite the Hamiltonian 
(after a Jordan- Wigner transformation) as 

JV N 

H =-2 E^- +1 + h.c] + e J2(-l) j (Sj + 1/2) 

.7=1 3=1 
N N 

= -tJ2lS*S* +1 + Spy +1 ] + e^(-l)^|, (5) 

3=1 3=1 

where as usual — S x ± iS y . The number operator rij 
in fermionic terminology is identified as Sj +1/2 in the 
spin terminology. 



B. Matrix Product States 

We refer to previous workEl : for a derivation of the 
ansatz and for details concerning the calculations in this 
section. A general matrix product state takes the form 



\Q)n = ^2tr[QA[s N ] ■ ■ ■ A[ Sl }]\s N ■ ■ ■ Sl ), 



(6) 



where Q is an m x m matrix containing the boundary 
conditions on the chain, A[s] is an m x m projection ma- 
trix obtained either from a DMRG calculation or varia- 
tionally, and Sj is the quantum number associated with 
site j. The projection matrix A contains the informa- 
tion about which states to keep when the lattice is aug- 
mented with one site. The number of degrees of freedom 
in A is reduced by preservation of orthonormal bases: 
J2 B A[s]A^[s] — l m . Further reduction of the number of 
free parameters is possible by exploiting symmetries of 
the system. In our variational calculations we have used 
particle-hole symmetry and conservation of the number 
of particles. Generally, the projection matrix A[s] is built 
up from states that form irreducible representations of 
the symmetry group of the Hamiltonian. 

In terms of the spin Hamiltonian Eq. (^), each factor 
A[s] in Eq. (g) adds a spin-1/2, hence taking a half- 
integer (hi) total spin into an integer (i) total spin and 
vice versa. This implies that we can define our projection 
matrix A with an off-diagonal block structure 



A[a] 



o Ai^iW 

Ai^Js] 



(7) 
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It is convenient to introduce the following mapping, 
denoted ~ from a local s x s matrix M to an m 2 x m 2 
matrix M : 



(8) 



Just using the block structure of A[s], one can showta 
that the eigenvalues of an operator M appear in pairs 
±A. We will frequently interpret the eigenvectors of M 
(of length m 2 ) as matrices of size m x m. 

The matrix 1, i.e. the ^-image of the identity ma- 
trix, plays an important role in the theory. Since 1 is 
guaranteed to have an eigenvalue 1, due to the orthonor- 
malization condition, there also exists an eigenvalue —1. 
The block structure of the projection matrix also implies 
that there may occur eigenvalues due to mixing of the 
integer and half-integer representations. However, these 
eigenvalues are spurious (unphysical) , in the sense that 
they do not affect the correlation functions, and can be 
removed simply by working with two different A matri- 
ces formed of the two off-diagonal blocks in Eq. (0). For 
this reason we leave these spurious eigenvalues aside in 
the subsequent discussion. 

For a translationally invariant system Eq. (^|) can be 
generalized to a Bloch-wave ansatz: 



|0, fc) 



N 



i,{s} 
x\s N 



ijk 



ti[A[s N }---A[s j+1 }Q---A[ Sl }] 



(9) 



where k is the momentum. 

The ground state of our model is translationally invari- 
ant. This implies that the matrix Q in Eq. (Q) should 
satisfy [Q, A[s]] — for all s, or_equivalently, Q must be 
a generalized right eigenvectorE3 of 1 with correspond- 
ing eigenvalue 1 (provided that this eigenvalue is non- 
degenerate). This means that Q ~ l m and hence our 
ground state ansatz takes the form, 



\1)n = ^2tx[A[s N ] ■ ■ ■ A[si]]\s N ■ ■ ■ s^. 



(10) 



Note thatpif we use a different sign convention in the 
HamiltonianE!3, the ground state would have momentum 
7r and hence we would have to choose a Q that anticom- 
mutes with A[s], {Q,A[s]} — for all s. From this it 
follows that we must choose Q as the generalized (right) 
eigenvector corresponding to the eigenvalue —1 of 1. Let 
us denote this generalized eigenvector by R for future 
purposes. This indicates that R is associated with mo- 
mentum 7r. Further evidence for this is given in con- 
nection to the discussion of the single particle excitation 
spectrum, see Section III C . 



C. Correlation functions 



Suppose we want to compute the ground state correla- 
tion function C(l) between two local operators Mj and 
acting on sites j and j + 1 respectively. Since we are 
working with a fermionic model it is necessary to distin- 
guish between local operators depending on whether they 
commute or anticommute on different sites. We refer 
to these operators as bosonic and fermionic respectively. 
For example, the density-density correlation is expressed 
in terms of two bosonic operators, while the particle- hole 
correlation is expressed in terms of two fermionic oper- 
ators. We will use the superscripts B and F to denote 
bosonic and fermionic operators respectively. For bosonic 
operators the correlation function is given by 



C(l) = (MjMj^i) = (l\l)- 1 tr[M B n l - 1 M B n 



N-l-ll 



(11) 



If we instead are interested in a correlation function be- 
tween fermionic operators, we have to keep track of the 
number of fermions between the sites j and j + 1 in order 
to get the phases correct. Defining F as a diagonal matrix 
with diagonal elements (—1, 1), we find the expression to 
be 

£7(0 = (M^M^) = (l\l)- 1 ix[^p- 1 M^»l If - 1 - 1 ]. 

(12) 



Eqs. (^lj) and ( |l2| ) implyQ that in general a correlati 
function takes the analytical form 



ion 



£7(0 = ^2a l sgn\ l i exp[-l/£ i ], 



(13) 



where £j = — 1 / In | A» | , and Aj are the eigenvalues of 1 or 
F, depending on the statistics of the operators. The a^s 
are coefficients that depend on the operators in the corre- 
lation function. Thus, the eigenvalues of Aj determine the 
possible correlation lengths in the system and it is there- 
fore important to investigate the spectrum of 1 and F. 
Eigenvalues that fulfill |A;| = 1 can potentially give rise 
to true long-range order. Due to normalization, 1 is guar- 
anteed to have eigenvalues ±1, which potentially could 
give long-range order in the bosonic correlation functions. 
In Appendix ^ we show that the spectrum of F differs 
from that of 1 only by a factor i. Hence, fermionic and 
bosonic operatorss have the same set of possible corre- 
lation lengths, which means that F has eigenvalues ±i 
and these can give rise to true long-range order in the 
fermionic correlation functions. Finally, we note that 
negative and imaginary eigenvalues correspond to oscil- 
lating correlation functions. 
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D. Overlap of DMRG states 

Suppose we perform two different DMRG calculations, 
keeping mi and m 2 states respectively. The overlap 
jv(lmi|lm 2 )iv ; where |l m3 )jv is the normalized ground 
state obtained by keeping 7772 states etc., can be com- 
puted as follows, 



w(lTOi|lm 2 )jv — tr[T^ l!m2 ], 



(14) 



where we have defined the mixed (miW2) X (7774771.2) ma- 
trix l mi as 



(15) 



In order to calculate the projection operators A[s], we 
have performed standard DMRG calculations by using a 
superblock of the form: 



the infinite lattice algorithm, and by adding a single site 
per iteration to each block. When we have an on-site 
potential e present in the problem we have to keep four 
projection matrices in order to completely describe the 
system, see Fig. [jj From these four matrices we form two 
projection operators A + and A~ taking us from positive 
to negative on-site potential and vice versa. Explicitly: 



.4 d 



A±\s] 



A?[a] 




Note that this overlap would be difficult to compute with- 
out the matrix product formalism since we have no map- 
ping between the different basis states of the two DMRG 
calculations due to renormalization. In contrast, all ma- 
trix product states are formulated in terms of the fixed 
{|sjv • • • Si)} basis (rather than renormalized basis sets), 
with the projection matrix just providing the amplitudes. 

Using Eq. (Ill) we find that the overlap decays expo- 
nentially as A^vwhere A is the (in absolute value) leading 
eigenvalue of l TOl , m2 . 

The overlap between DMRG states for different num- 
ber of kept states gives a measure of the gain in accuracy 
obtained by increasing the number of kept states. We 
consider this measure to be more universal and relevant 
than the usual measure used, namely the truncation of 
the density matrix, 1 — trp e , where p e is the truncated 
density matrix. The old measure is algorithm-dependent 
and can, for instance, be made equal to zerca, despite the 
fact that an exact calculation is not p erformed, simply 
by using a superblock configuration | b j 





FIG. 1. The figure shows the four projection matrices 
needed to describe the system. Filled discs denote sites with 
positive on-site potential and circles denote sites with nega- 
tive on-site potential, hi denotes half-integer representations 
and i denotes integer representations. Andersson et. al. 

In each DMRG iteration we update either the hi — * i 
or the i — > hi matrices. 



E. Computational methods 

Due to the large dimensions of l TOlir7l2! namely 
(7711777,2) x (7711777,2)7 it becomes necessary to use itera- 
tive eigenvalue routines that require no explicit construc- 
tion or storage of l mi , m2 . Moreover, l mijTri2 iSpCion- 
symmetric. We have used the Arnoldi algorithmEd to 
handle these problems. 

Furthermore, the computations become much more ef- 
ficient if we rewrite the operation of l mi , m2 on an (77117712) 
vector v as a matrix product with 71 interpreted as an 
mi x TU2 matrix: 



In this way we only need to operate with mi x 7772 matri- 
ces, and the eigenvalues of l mi , m2 can easily be obtained. 



III. RESULTS AND DISCUSSION 
A. Convergence of the DMRG 

In this section we will discuss the convergence of the 
DMRG. First, we demonstrate that the DMRG projec- 
tion operator of the critical system converges to a fixed- 
point, and thus, justifying the matrix product ansatz 
when studying the thermodynamic limit of the DMRG. 
Second, we check the convergence of the ground state 
with respect to the number of kept states by using the 
overlap measure. 

The fundamental assumption of the matrix product 
approach is that the projection matrix converges to a 
fixed-point with respect to N, i.e. limjv— >oo Atv[s] = A[s]. 
In order to show this, we define the matrix norm || • \\ max 
via, 



\A\ 



max — max 



\A 



(16) 
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In addition, we make a consistent enumeration (with re- 
spect to quantum numbers) and use a fixed sign conven- 
tion of the states in the system blocks. It is then easy 
to study the convergence of the projection operator by 
measuring the quantity r(N) = (||Ajv+i — Ajvilmoa; + 
— Aff_i\\ max )/2, where N is the number of DMRG 
iterations. In Fig. || we have shown results from such 
calculations. From the figure we see that the conver- 
gence of the projection operator seems to be exponential 
with respect to the number of DMRG iterations and that 
the convergence rate decreases when the number of kept 
states is increased. We have also found that the ground 
state energy per site converges much faster than the pro- 
jection matrices and is therefore not a good indicator on 
whether or not a fixed-point has been reached. 











(36,32) 


\ 

\(4,6) 

\ 


(12,12) 


(21,20) 



50 100 150 200 



DMRG iterations, N 
FIG. 2. The norm r(N) is shown as a function of the num- 
ber of DMRG iterations, N, for the gapless case. It is clear 
from the figure that the DMRG projection operator converges 
with respect to N. Furthermore, the convergence seems to be 
exponential with respect to N. The peaks in the (36, 32) curve 
indicate that either the DMRG has changed the states kept 
in the Hilbert space basis or our sign-fixing procedure of the 
states has failed. The notation (7711,7712) means mi states in 
the integer representation and 7712 states in the half-integer 
representation. 



(mi,m 2 ) (4,6) (12,12) (21,20) (36,32) 
(4,6) 0.00231 0.00328 0.00381 
(12,12) 0.00231 0.000277 0.000644 
(21,20) 0.00328 0.000277 0.000134 
(36,32) 0.00381 0.000644 0.000134 

TABLE I. Leading eigenvalues 1 — A, govering the overlap 
between different DMRG ground state wavefunctions. The 
notation (mi, 7712) means mi states in the integer representa- 
tion and 7712 states in the half-integer representation. 

In Table. I we have presented the eigenvalues that 
govern the overlap in the thermodynamic limit between 
ground states obtained by keeping different numbers of 



states. We have chosen to write out 1 — A instead of A 
since this gives a more direct measure of the error. 



B. The spectrum of 1 and correlation lengths 

Since the DMRG projection operator converges to a 
fixed-point, the correlation functions of the DMRG in 
the thermodynamic limit are given by Eq. (£§) and the 
correlation lengths are determined by the eigenvalues of 
1 (the correlation lengths obtained from F are identical) . 
An analysis of the spectrum of 1 is therefore pivotal. 

We have found that the eigenvalues ±1 are non- 
degenerate, and that all the other eigenvalues fulfill 
A I < 1. Only the eigenvalues ±1 can give rise to infi- 
nite correlations lengths. However, it turns out that that 
the density-density and the particle-hole operators are 
orthogonal to the corresponding eigenvectors, and hence 
the correlation lengths will be determined by other eigen- 
values. 

An interesting question is which local operators M that 
potentially can give true long-range order in the correla- 
tion functions Eqs. ([□]) and jl^). We show in Appendix 
[b| that true long-range order for bosonic operators is not 
possible if tr[M B ] = 0. For fermionic operators we find 
that true long-range order is not possible for off-diagonal 
operators. This explains why there is no true long-range 
order in the density-density and particle-hole correlation 
functions. The proof exploits that the Hamiltonian con- 
serves the number of particles and that it is particle- 
hole symmetric. If we break the particle-hole symmetry, 
we are only guaranteed that the particle-hole correlation 
function can not be truly long-ranged. 

Thus, the DMRG will approximate infinite correlation 
lengths by finite, and it is interesting to investigate how 
the DMRG correlation lengths depend on the number of 
kept states and on the gap of the system. 

First of all we need to identify the leading eigenvalues 
of F and 1 govering the particle-hole and density-density 
correlations respectively. These eigenvalues can be iden- 
tified either by a matrix product calculation of the respec- 
tive correlation function, or by measuring the correlation 
length directly in the DMRG calculation. Degeneracies 
in the spectrum are also instrumental in identifying the 
leading eigenvalue. For instance, if we want to compute 
the particle- hole correlation function, we expect this cor- 
relation to couple to an eigenvalue that is two-fold de- 
generate since the hole-particle correlation function has 
an equal correlation length. The density-density correla- 
tion function, on the other hand, will couple to a non- 
degenerate eigenvalue since there is no symmetry related 
correlation function that demand an equal eigenvalue. 

Using the analytical result for the correlation length 
given by Eq. (||) together with A = exp[— we find 
that the exact expression for the eigenvalue dominating 
the particle-hole correlation is: 
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^(M) = ^ + J--. (17) 

The * is used to indicate that this is an exact value. The 
expression will be used as a reference when we evaluate 
our numerical data. 




0.85 1 1 1 1 1 ' 1 ' 1 1 1 

0.00 0.02 0.04 0.06 0.08 0.10 

E 

FIG. 3. The eigenvalues of F govering the parti- 
cle-hole correlation length versus e for different numbers of 
kept states. The solid line corresponds to the exact re- 
sult from Eq. ( |l7| ) From the bottom to the top, the point 
sets correspond to the following numbers of kept states: 
(4, 6), (12, 12), (21, 20), (36, 32). Due to numerical problems, 
we can not continue the set with mi or m,2 > 12 to larger e 
values. 

10° t , , 



number of kept states in our truncated Hilbert space, 
the accuracy of the correlation length increases. For the 
case e = we see that the eigenvalue \X p h\ approaches 
the exact value 1, i.e. an infinite correlation length, as 
the number of states is increased. The convergence of 
\\ p h\ towards 1 is more clearly seen in Fig. 0, where we 
consider the gapless case and plot 1 — | A p ^" versus the 
number of kept states, m = mi + m-i . Thus, we conclude 
that the DMRG gives exponentially decaying correlation 
functions in the gapless case, but as the number of states 
is increased the correlation length grows towards infinity. 

We can actually make this conclusion more quantita- 
tive. As is seen in Fig. ||, the eigenvalue \\ v h\ behaves as 
\Xph\ — 1 — km~" . Thus the correlation length behaves 
as 

* - mil- 1 *™-*! * l mP - (18) 

That is, the correlation length scales as a power of m. We 
find the exponent (3 ~ 1.3 and k ~ 0.45. The density- 
density correlation function gives similar results. 

In any realistic DMRG computation, the correlation 
function will be given by a sum of a finite (even if large) 
number of exponentially decaying functions according to 
Eq. (|l5|). Keeping only a few states, we have seen that 
the correlation function approximates the true power-law 
(Eq. (^)) for short correlations, but as we increase / we 
eventually end up with an exponential decay. Increasing 
the number of kept states will make the correlation func- 
tion look like a power-law for rather large I, but in the 
end, when I — > oo, it will always behave as an exponen- 
tially decaying function with correlation length given by 
Eq. ©. ' 



C. The single-particle excitation spectrum 
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FIG. 4. Convergence of the eigenvalue A of F govering the 
particle-hole correlation length. In the figure we show 1 — |A| 
versus the number of kept states, m = mi +rri2, in the gapless 
(e = 0) case. 

In Fig. U our results for the particle-hole correlation 
length are shown. It is clear that as we increase the 



In order to study the single particle excitation spec- 
trum, we have used the Bloch-wave ansatz of Eq. (||) 
and the pole-expansion techniques to calculate the spec- 
trum. The result is shown in Fig. |[ The curve shows a 
pair of excitations corresponding to a single particle/hole 
(or spin S z = ±1). We can see that the dispersion re- 
lation obtained from the Bloch-wave ansatz is in good 
agreement with the exact dispersion E(k) = sin k, ex- 
cept close to the Fermi points, where the gap closes. In- 
stead of having a linear form, the calculated dispersion 
relation has the form E(k) ~ Ao + v^k 2 close to k = 0. 
Furthermore, our excitations have negative energies close 
to the Fermi points, a consequence of a defect gijcomd 
state. Such a negative energy gap was also foundij to 
appear for the biquadratic spin-1 chain somewhere be- 
tween the Heisenberg point and the Takhtajan-Babujian 
point. We have investigated how the size of this nega- 
tive energy gap depends on the number of kept states in 
the truncated Hilbert space. These calculations are com- 
putationally demanding and also sensitive to numerical 
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errors. However, the results we have indicate that the 
size of the negative energy gap decreases as the number 
of states is increased, but the numerical control was poor 
for mi + m 2 > 12. 

A possible explanation of this defect ground state could 
be that the DMRG has an instability of some sort. Wc 
have investigated whether the DMRG at e = is unstable 
against breaking the translational symmetry of ground 
state by starting the DMRG calculation with a non-zero 
staggered on-site potential and then, after about 20 iter- 
ations, turning the potential off. We then let the projec- 
tion operator converge and compare it with the projec- 
tion operator obtained by a DMRG calculation with the 
on-site potential turned off all the time. However, the two 
converge to the same limit (within the numerical accu- 
racy) with only one exception. If we keep 3 states in the 
integer representation and 6 states in the half-integer rep- 
resentation, we find that the DMRG is actually unstable 
and we obtain an energetically more favorable state by 
breaking the translational symmetry. However, it is suffi- 
cient to add a single state in the integer representation in 
order to remove this instability. In addition, we have also 
performed variational calculations allowing for a ground 
state with periodicity two, but the energetically lowest 
state turns out to be translationally invariant. Thus, it 
seems like the DMRG is stable against forming a ground 
state which is not translationally invariant. 
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Bloch wave 
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FIG. 5. Single particle energy dispersion relation E(k). 
The Bloch-wave result deviates from the exact result only 
close to the Fermi points, (mi, 7712) = (4,4). 

There is an interesting symmetry concerning the dis- 
persion relation in Fig. Q We find that the dispersion 
relation has the symmetry E(it — k) = E{k) in the ther- 
modynamic limit. This is in fact a consequence of the 
block structure of A[s], and should therefore be a char- 
acteristic feature of many half-integer spin systems. To 
prove this we use the R matrix to explicitly construct a 
Bloch state of momentum k + 7r from a state of momen- 
tum k and show that these two states have equal energy. 



Thus, we have shown that E(k) = E(k+Tr), but this also 
proves that E{k) = E(ir — k) since we have the sequence 
of mappings: 

E(k) E(-k) E(tt - k), (19) 

where V is a parity transformation, an exact symmetry 
of our model. The details of the proof can be found in 
Appendix |c[ This symmetry is true in general for trans- 
lationally invariant half-integer spin systems with local 
interactions. As a test, we have checked numerically that 
it is true also for the isotropic antiferromagnetic spin-1/2 
Heisenberg chain. Furthermore, calculations on the spin- 
1 chain, as expected, lack this symmetry. The symme- 
try has recently been studied on more general grounds 
by Kladkdij, without any reference to the Bloch-wave 
ansatz. Similarly to our proof for the Bloch-wave ansatz, 
Kladko explicitly constructs a state with momentum fc+7r 
from a state with momentum k and then shows that these 
states have equal energy. 

IV. CONCLUSIONS 

Wc have investigated fundamental properties of the 
DMRG when applied to a gapless system of free fermions. 
We find that the DMRG projection operator converges 
to a fixed-point. This convergence means that states 
of a matrix-product form are identical to the DMRG 
states in the thermodynamic limit. By using the matrix- 
product formalism, we have found that DMRG calcula- 
tions give qualitatively wrong particle-hole and density- 
density correlation functions: the DMRG correlations de- 
cay exponentially, while the true correlations decay alge- 
braically. However, for short distances, the DMRG cor- 
relation function agrees with the exact result. The finite 
correlation length of the DMRG particle-hole correlation 
function scales as £ ~ m 13 , where m is the number of 
kept states. In addition, we have derived conditions for 
whether a general operator potentially can give rise to 
truly long range correlations or not. These conditions 
are found to be determined by symmetries of the Hamil- 
tonian. 

We have demonstrated that the matrix product formal- 
ism can be used to calculate overlaps between differently 
renormalized states. This makes it possible to directly 
compare states obtained from DMRG calculations, using 
different numbers of kept states. We propose this overlap 
as a criterion of convergence of DMRG states. 

From the matrix-product ansatz we obtain accurate 
values for the ground state energy. Furthermore, using a 
Bloch-wave ansatz, we find a dispersion relation for the 
excitations that is close to the exact result. Despite this, 
close to the Fermi points, where the gap closes, our exci- 
tations have negative energies, indicating that the ground 
state is defect. We have not yet been able to trace the 
origin of these negative energy excitations, although our 
calculations indicate that the magnitude of this negative 
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energy gap decreases as the number of kept states is in- 
creased. In addition, we have shown that the Bloch-wave 
ansatz for the excitation spectrum exhibits the symmetry 
E(k) = E(ir — k) for translationally invariant half-integer 
spin systems with local interactions. 
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APPENDIX A: SPECTRUM OF F 

Since F determines the possible correlation lengths of 
fcrmionic operators, it is important to understand the 
eigenvalue spectrum of this operator. We will in this 
appendix show that the eigenvalues of F are related to 
those of 1 by a factor i. That is, Xf = iAj.. To show this, 
we will start by constructing an eigenvector of F with 
eigenvalue i. This-construction is similar to the one used 
by Roman et. al£n. Using F[s,s'] — 5 s , s ii 2s+1 , we may 
write the matrix elements of the operator F as 

p(Y 2 ,m' 2 )(Y 1 ,m' 1 ),(-f2,m 2 )(-(i,m 1 ) 

— j2(m2-m 2 ) + l^(72 ! m 2)(72,"i2)j s j^(7i,"ii)(7l,™l)j s j ^ 

s 

where we have used that since the projection op- 
erator conserves the particle number, the element 
Ah ,m ),(7,«i)js] i s zer0 UI1 less m' = m + s. Next, we 
define the vector \ UF )^ 2 ' m2 ^ mi) = 5 l2ill 5 m2 . mi i 2m \ 
We will now show that \u F ) is an eigenvector of F with 
eigenvalue i. We have 



(F\u F )) ( ' 1 ' 2 ' m ' 2){l ' 1 ' m ' l) 

— j2(m 2 -m 2 )+l^(72> m 2)(72,m 3 )f s 



E 

71,72 mi, rri2,s 

xA(<< m '^< m ^[s}i 2m *5 mum2 5^ 2 



^4(72 :™2 ) (71 ,mi ) M ,™i ) (71 ,mi ) J 



VI 



2ml 



71 ,mi ,s 



which proves the claim. In the third line we have used 
the identity ^ s ^4[s]A T [s] = 1. Let us now show that the 

entire eigenvalue spectrum of F is related to that of 1 by 
a factor i. First of all we note that | detu^l = 1, which 
implies that the inverse Up 1 exists. In fact Up 1 = u F , i.e. 
uf is unitary. Furthermore, Uf satisfies the equation 

iu F A[s] = F[s, s]A[s]u F , (Al) 



since this equation is equivalent to F\up) = i\up), as 
can be seen by simply multiplying Eq. ( [Al] ) by [s] and 
summing over s. Define the unitary operator £ = (ttF<S>l) 
and consider 



^l£ f = i{uF ® 1) ( E ® A i s ] j ( U F ® 1) 

^ S ' 

= i^2(u F A[s]u F ) ®A{s] 

S 

= F[s, s]A[s] ® A[s](u F u F ®l)=F. 



This implies that F has exactly the same spectrum as il, 
whic h was our claim. In the third line we have used Eq. 



APPENDIX B: CONDITIONS FOR TRUE 
LONG-RANGE ORDER 

In this appendix, we show how symmetries of the 
Hamiltonian determine which local operators that poten- 
tially can give true long-range order in correlation func- 
tions. In order for a local bosonic operator M B to give 
true long-range order it must hold that at least one of 
the following expectation values is non-zero: (1|M S |1), 
(1\M^\R) and (R\M^\R), where |1) and \R) denote the 
eigenvectors of 1 with eigenvalues 1 and —1 respectively. 
These expectation values will be determined if we can 
determine the expectation values of A*[s] £g> A[s'] for all 
combinations of s and s' . We will frequently interpret the 
m 2 state- vectors |1) and \R) as m x m matrices denoted 
by 1 and R. The matrix R has the block- form 1 © — 1 
which implies that R 2 = 1. A subscript L on the matrix 
denotes that it represents the left eigenvector. Further- 
more, we will consider the projection matrix A[s] to be 
real. 

Let us start by relating different expectation values. 
We have 



T V1 



{R\A[s] ® A[s']\R) = tr[RlA[s]RA 

= -tr[RiRA[s]A T [s'}} = -tr[llA[s]A T [s']} 
= -{l\A[a]®A[ a r\\l), 

and similarly, 

(R\A[s] ® A[s']\l) = -{l\A[s] ® A[s']\R). 

The indices of A[s) correspond to different states and 
we will label these states as I7, m), where m is the particle 
number measured from half-filling and 7 is an integer la- 
beling particle-hole representations. The transformation 
of the state I7, m) under a particle- hole transformation 
B is 



B\j,m) = <f)('y,m)\j,-m) 



(Bl) 
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where 0(7, m) = ±1. That is, the state |7, m) trans- 
forms within the 7 representation of the particle-hole 
symmetry group. Note also that B 2 = 1, which im- 
plies that <j) is independent of to, since B 2 |7,to) = 
0(7, m)<£(7, -to)|7,to). 



We are now going to show that 1 



(7 / ,m'),(7,m) 



is zero un- 
less to' = to. To see this we write the equation (1| = (1|1 

as 



1 



(71, mi), (72, m 2 ) 



^(73,m 3 ), (74,77x4) 

s 73,74,7713,7714 
x ^(73,m 3 ), (71, mi) |" g j^(74, m 4 ), (72, rn 2 ) J s j 



Using the fact that A conserves the particle num- 
ber (see Appendix |X|), we conclude that mi — TO2 = 
TO3 — TO4. Since the difference between the to- values of 
is conserved under the action of 1, we 



a 



(71, mi), (72, m 2 ) 



can write \l = J-x, + where a^ only has non-zero 
matrix elements between states of equal particle num- 
ber and V L est is the remainder, mixing particle num- 
bers. This decomposition will not be mixed under the 
action of 1 and hence both matrices must be eigenstates 
of 1 with eigenvalue 1. This eigenvalue is, however, non- 
degenerate, which means that either of 1 L and V L est is 
zero. Since tra^ = (1|1) = 1, must contain t L and 
hence V L est = 0. This means that 1l only connects states 
containing an equal number of particles. 

Furthermore, since Rl = iz,-R and R also conserves 
the number of particles (it is diagonal), we conclude that 
Rl conserves the number of particles. 

Let us now consider the following expectation value: 

(l\A[s] ® A[s']\l) = tr[ll A[s}A T [s'}} 

— ^^(71, mi), (72,7712) 

{li,mi} 

x J 4(72,m 2 ),(73,m 3 ) ^^(71,7711), (73,7713) ryj 

Using particle number conservation of A and a l we con- 
clude that the expectation value is zero unless s = s', 



(l\A[s]®A[s']\l)~6 s 



(B2) 



Similarly, using Rl instead of a^, one finds 

(l\A[s]®A[s'}\R) ~6 S , S ,. (B3) 



To summarize, Eqs. (B2) and flB3|) follow from the 
particle-number conserving property of the Hamiltonian 
together with the uniqueness of the eigenvalue 1 of 1. 
In the next paragraph, we will investigate the expecta- 
tion values when s = s'. Note that the Pauli matrix 03 
and a form a complete basis for all diagonal operators 
M B , which according to the above results are the only 
operators that can give true long-range order. 

Recall the defining relation of the projection operator: 



W,m')= ^ (7 ' VM7 ' m) [^]|7,m)®| Sj ). (B4) 

(7,m),Sj 

Applying the particle-hole transformation to the defin- 
ing relation Eq. (B4), we find that A has the following 
symmetry: 

A (7',m'),(7,7n)[ Sj .] = 0( 7 ')0( 7 ) (_ 1)^(7' ,-m'),(7,-m) r_^ _ 

(B5) 

Let us write down the operator form of A[s] 

A[ Sj )= Y, ^ (7 '* m ' )l(7,m) fe]|7',m')(7,m|. 

(7', m'), (7,771) 

Applying a particle-hole t ran sformation to A[s], making 
use of the symmetry Eq. (BE), one finds that 



BA[ Sj ]B = (-l) j A[- Sj 



(B6) 



Using the definition of 1 we find the transformation prop- 
erty of 1 under particle- hole transformations: 

(B ® B) Y A[s] ® A[s] (B ® B) = ^ A[-s] ® A[-s] . 



That is, 1 is invariant under particle-hole transforma- 
tions. Let us also consider the transformation of 03, 
where 173 is a Pauli matrix. Using exactly the same tech- 
nique, we find that (B ® B)ai{B ® B) — —03. 

(1| and (1\(B ® B) must both have eigenvalue 1 of 1, 
since 1 is invariant under particle-hole transformations. 
Using the non-degeneracy of this eigenvalue, we conclude 
that (1\(B®B) = e- l9 (l\ and similarly (B<g>B)|l) = e l9 \l). 

Let us now compute the expectation value of 03, 

(1|SJ|1) = (1|(B ® B) 2 ^ 3 (B ® B) 2 \l) = -(1|S5|1), 
and thus (ljirjjll) = 0. Since (1|1|1) = 1 we arrive at 

(l\A[s] ® A[s']\l) = -(R\A[ s }®A[s']\R) = ^y. ( B7 ) 

It only remains to consider expectation values of 1 and 
between (1| and \R), Trivially, (1\1\R) = 0, since (1| 
and \R) are eigenvectors of 1 with different eigenvalues. 
Using the structure of |1) and \R), let us show that they 
must transform with the same phase factor. Noting that 
B must be block-diagonal, it follows that [B, R] = 0. As- 
suming (B <8>B)|1> = BB T = e' e \l) we find 

(B ® B)\R) = BRB T = RBB T = Re t9 t = e l6 \R). 
Using this property we find 

(\\ai\R) = -(1\*S\R)=0 
and we conclude 

(l\A[a] <8> A[s']\R) = -{R\A[s] <8> A[s']\l) = 0. (B8) 
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Using the derived expressions for the expectation val- 
ues, a general bosonic operator M B can couple to the 
eigenstates with corresponding eigenvalues ±1 only if 
trM B ^ 0. In particular, the density-density correlation 
function can not be truly long-ranged when the Hamilto- 
nian has particle-hole symmetry and conserves the num- 
ber of particles. 

Let us also investigate correlation functions between 
fermionic operators, M F , in order to be able to draw con- 
clusions concerning the particle-hole correlation function. 
We will only show that off-diagonal fermionic operators 
can not give rise to truly long-ranged correlation func- 
tions. From Eq. (|l2|) we see that true long-range order 
is possible only if some expectation value of the form 
is non-zero, that is S + must connect eigen- 
states of 1 and F with the corresponding eigenvalues 
having absolute value 1. We will now show that this is 
impossible due to the particle-number conserving prop- 
erty of A[s]. Note that all the important eigenvectors are 
proportional to S m , m '- Hence, an expectation value of 
the above form will be 

(l\S + \up) ~ ^""^ Sm,m' &m, 1 /2+m" S m ' , - 1/2+m'" $m" ,m"' 

= o, 

where the second and third delta-functions come from 
the particle- number conserving property of A[s]. Thus 
we can not have true long-range order in the particle- 
hole correlation function. 

Note that we have not proved that the correlation func- 
tion between two traceless fermionic operators (like 03) 
can not be long ranged. On the contrary, this is the 
structure of the string order correlation function in the 
spin-1 chain, which is long-ranged. 

APPENDIX C: PROOF OF THE BLOCH-STATE 
SYMMETRY 

In this appendix, we prove that the energy spectrum 
obtained from the Bloch-state ansatz exhibits the sym- 
metry E(k) = E{tt — k). Since the symmetry E(k) 
I follows from parity being a good quantum num- 

ber, we only need to show that E(k) — E(k + tt). The 
strategy used in the proof is to construct a Bloch-state 
of momentum k + tt from a state of momentum k and 
to show that expectation values of these two states are 
equal in the thermodynamic limit. 

Let us first recall some properties of the eigenvector R 
of 1. We will use the convention that \R) is an m 2 vector 
and R is an m x m matrix and similarly we will write the 
eigenvector of 1 corresponding to the eigenvalue 1 as |1) 
or 1. Using the block-diagonal structure of R it follows 
that {(1 <8> R), M} — and that [(R (g> R), M] = for all 
local operators M. 

We also need the property (R <g> R) lim^oo 1™ = 
linin^oo 1™. To show this we recall that all but two eigen- 



values of 1 have absolute value less than one. The corre- 
sponding eigenvectors will be annihilated by linin^oo 1™. 
Hence we may write 

lim T" = |l)(l| + (-l)"|i?)(i?|. 

n — >oo 

If we let this operator act on a general state \ip) = 
J2i we obtain 

lim l n \^) =Tj h \l) + (-l) n i> R \R). 

n — >oo 

Now, if we act on the resulting vector with (R ® R) , we 
obtain (now we use the matrix form of the vectors) 

(R®R) lim 1*|V>) = (R<8R)[1>i\l) + (-l) n ip R \R)] 

n — >-oo 

= i) 1 RtR T + {-l) n ^ R RRR T 

= i/>il + {-l) n ip R R = lim T n |V) • 

n — >oo 

Thus we have shown that (R eg) R) lim^-joo 1™ and 
linin^oo 1" act equally on a general state and hence the 
operators must be identical. 

We are now ready to show that for all local operators 
M, it holds that 

LHS = (RQ',k\M\RQ,k) 

= (Q', k + n\M\Q, k + n) = RHS , (CI) 

where the Bloch-states are defined in Eq. (||) . If we apply 
this result to the Hamiltonian operator and the normal- 
ization, we obtain the result that the states \RQ, k) and 
\Q,k + tt) have equal energy. That is, we have a map- 
ping from a state of momentum Ho a state of momen- 
tum k + n with equal energy, which proves the symmetry 
E(k) = E(k + tt). 

In order to prove Eq. ( |ClD , we begin by writing the 
Bloch-state \RQ,k) as 

\RQ,k)= e lj(k+7T) tr[RA[ SN }---A[ S3+1 }Q---Al Sl }] 

J.fsj} 

x|stv • • ■ Si) , 

where we have used {i?,yl[s]} = to mcpe the R to the 
left side of t he t race. Thus we may writdlJ the left-hand 
side of Eq. (|Cl]) as 

LHS = e y ' (fc+7r) e-' y "( fe+7r ) x tr[(R g> R)l N ~ j 
ho' 

x(U ® Qp^Ml^'^iQ'* ® 1)V] , 

where we have assumed that M acts on the site I. Now, 
when we go to the thermodynamic limit N — > 00, we will 
always have a factor I 00 somewhere in the trace. If we 
could move (R (g> R) through the trace until it reaches 
the I 00 factor, the (R <g> R) would be annihilated and we 
would be left with an expression that is 
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eiKfc+^e-iAfc+^trp^-^l ® Q) 
xP'-'mT'^'-^Q'* ® 1)D''] = RHS 

and our proof would be complete. It turns out that it 
is always possible to perform such a move. If the factor 
1°° is not between the factors (1 ® Q) and (Q'* ® 1), 
we just commute (i? ® R) through the ^-operators until 
it reaches 1°° and gets annihilated. If the factor 1°° is 
between (1 <8> Q) and (Q'* ® II) we can make the split 
(R <g> R) = (1 <8> R)(R <g> 1) and move these factors in 
different directions until they meet between (1 <g> Q) and 
(Q 1 * Cg) 1) and get annihilated by 1°°. During the move- 
ments we will pick up a factor (— 1) J +( N ~i ), and since 
we are only considering chains of even length, this factor 
is equal to 1. Note that the proof works even if the oper- 
ator M does not act on a single site, but, as for example 
the Hamiltonian, acts on a couple of neighboring sites. 
The important thing is that the operator is local, so that 
a factor of 1°° always can be found in the trace. 
This completes the proof. 
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